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The bifurcation and chaotic behaviour of unidirectionally coupled deterministic ratchets is studied 
as a function of the driving force amplitude (a) and frequency (oj). A classification of the various 
types of bifurcations likely to be encountered in this system was done by examining the stability of 
the steady state in linear response as well as constructing a two-parameter phase diagram in the (a — 
uj) plane. Numerical explorations revealed varieties of bifurcation sequences including quasiperiodic 
route to chaos. Besides, the familiar period-doubling and crises route to chaos exhibited by the one- 
dimensional ratchet were also found. In addition, the coupled ratchets display symmetry-breaking, 
saddle-nodes and bubbles of bifurcations. Chaotic behaviour is characterized by using the sensitivity 
to initial condition as well as the Lyapunov exponent spectrum; while a perusal of the phase space 
projected in the Poincare cross-section confirms some of the striking features. 
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I. INTRODUCTION 

Coupled nonlinear systems possess a rich catalog of 
potentially useful dynamical behaviours including bifur- 
cations, chaos and synchronization. They are central to 
the understanding of a wide variety of extended systems, 
e.g. a line of lattice oscillations or coupled plasma wave 
modes. A system of two oscillators could make a Hopf 
bifurcation to a second incommensurate frequency and 
follow a quasiperiodic route to chaos, in addition to pe- 
riod doubling and intermittency. Such routes to chaos 
was first observed by Ruelle and TakenspJ. Synchroniza- 
tion phenomena in coupled or driven nonlinear oscillators 
are of fundamental importance and have been extensively 
investigated both theoretically and experimentally in the 
context of many specific problems arising in laser dynam- 
ics, electronic circuits, secure communications and time 
series analysis |2( to mention a few. In a recent study, we 
observed phase synchronization in both unidirectionally 
@ and bidirectionally 3 coupled deterministic ratchets 
that exhibits intermittent chaos. In it was shown 

that the transition to the synchronous regime is charac- 
terized by an interior crises transition of the attractor in 
the phase space. Specifically two analytic tests based 
on the Fujisaka and Yamada approach and that of 

Gauthier and Bienfang were employed to verify the 
stability of the synchronous state in reference 0]. 

A clear understanding of coupled oscillator systems 
requires a study of the full spectrum of the operating 
regimes, including the nonsynchronous state |||. It has 
been reported that a group of interacting chaotic units 
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that exhibits synchronization show some bifurcation cas- 
cades from disorder to partial and global order when 
varying the coupling strength Ding and Yang [To| 
showed the existence of intermingled basins in coupled 
Duffing oscillators that exhibit synchronized chaos and 
conjectured that intermingled basins can be easily real- 
ized in the context of coupled oscillators and synchro- 
nized chaos. Extensive bifurcation analysis of two cou- 
pled periodically driven Duffing oscillators was done by 
Kozlowski et al. They showed that the global pat- 

tern of bifurcation curves in parameter space consists of 
repeated subpatterns similar to the superstructure ob- 
served for single, periodically driven, strictly dissipative 
oscillators. This study was recently extended to two cou- 
pled periodically driven double-well Duffing oscillators 
by Kenfack 01 • The results revealed a striking depar- 
ture from the single-well Duffing oscillators studied by 
Kozlowski et al [ill ]. 

The literature is abundant with studies on coupled os- 
cillator systems. At this point, we wish to refer the reader 
to some other recent studies that are related to the sub- 
ject of this work [13j]-[2(|- Motivated by this series of 
studies, we aim in this paper to explore some dynamical 
features hitherto not fully explored in the coupled ratch- 
ets with emphasies on the bifurcation structures proceed- 
ing the stable synchronous regime. The rest of the paper 
is organized as follows: Section 2 describes the equations 
of motion of the chaotic ratchet model and comment on 
the synchronization behaviour of the coupled system that 
is of interest in this study. In Section 3, we carry out a 
linear stability analysis of the system projected on the 
Poincare map, while in Section 4, numerical results of 
bifurcations are presented. Chaotic behaviour is charac- 
terized in Section 5. The paper is concluded in Section 
6. 
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II. THE CHAOTIC RATCHET MODEL 

Let us consider the one-dimensional problem of a par- 
ticle driven by a periodic time-dependent external force 
under the influence of an asymmetric potential of the 
ratchet type [2l]]-|2!|. The time average of the external 
force is zero. In the absence of stochastic noise, the dy- 
namics is exclusively deterministic. The dimensionless 
equation of motion for a particle of unit mass moving in 
the ratchet potential V(x) is given by 



bx 



dV(x) 
dx 



a cos(woi) 



(1) 



Where time t has been normalized in the unit of u>q , 
the small resonant frequency of the system. Here, a and 
b are the forcing strength and the damping parameter, 
respectively. The dimensionless potential V{x) is given 
by 



V(x) = C- 



1 



4tt 2 S [ 



i27r(a;-xo)+0.25sin47r(a;-xo)], (2) 



with constants C ~ 0.0173 and 5 ~ 1.6. The potential 
V{x) is shifted by a value xq in order to place its mini- 
mum at the origin (see fig. Notice that this asymmet- 
ric potential is periodic and has infinite potential wells. 

The extended phase space in which the dynamic is tak- 
ing place is three-dimensional, since we are dealing with 
an inhomogeneous differential equation with an explicit 
time dependence. We can rewrite equation (1) in au- 
tonomous form as a three dimensional dynamical system 
described by the coordinates: 



y, 

a cos z — by 



dV(x) 
dx 



(3) 



Z = U! D . 



Since equation (3) is nonlinear its solutions allow the pos- 
sibility of periodic and chaotic orbits. For the purpose of 
this study, we consider a system of two identical unidi- 
rectionally coupled chaotic ratchets governed by: 



xi + bx\ 



dV(x! 
dx\ 



a cos(ujt) 



(4) 



dV{x2) 

X2 + 0x2 H ; = a cos{ujt) + c[x\ — X2) (5) 

dX2 

Where c is the coupling parameter. For a — 0.08092844, 
b = 0.1 and u> = 0.67, systems (4) and (5) exhibit sta- 
ble phase synchronization when 0.89 < c < 1.05. The 
synchronization dynamics for c = 0.95 shown in fig. [3 is 
determined by the behaviour of the following difference 
x = X2 — x\. The phase difference is locked (full syn- 
chronization, x = 0) after a short transient time. Prior 
to the phase locking event, the strange attractor under- 
goes a crises transition during which it is gradually de- 
stroyed and as the synchronization regime is approached, 



the attractor is again re-built Q ■ To obtain approximate 
visualization of the attractors and their bifurcations we 
follow the procedure described in 0,0]. That is, we in- 
vestigate the dynamics in the Poincare sections denoted 

byE- 




FIG. 1: The ratchet potential (C = 0.0173, S = 1.6) 



III. STABILITY ANALYSIS 

To begin with, let us consider the basic mathematical 
principle underlying some expected topological structure 
of the coupled ratchets. Here we carry out a stability 
analysis of the coupled system (4) and (5) which can be 
written equivalently as 



(6) 



Where V{x\,Vi,X2,V2,u) is an autonomous vector field 
and ip(a,b,c,u>) is an element of the parameter space. 
Thus (4) and (5) become: 



1000 1500 2000 2500 3000 

time 



FIG. 2: Synchronization dynamics showing the quantity 
x = X2 — xi as function of time, for c = 0.95, b — 0.1, 
a = 0.08092844 and u = 0.67 
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xi = Vi 



v\ = a cos(ujt) — bv\ H 

ind 

X2 = v 2 



i>2 = a cos(ujt) — bv2 + 



2 cos 2ir(xi — xq) + cos 4ir(xi — xq) 

2 cos27r(a;2 — xq) + cos47r(a;2 — xq) 
I 



+ c(vi - v 2 ) 



(7) 



The system described by (7) generates a now cf) = {4> T } 
on the phase space K 4 x S 11 such that a global map 



exist, with 9q being a constant determining the lo- 
cation of the Poincare cross-section defined by — 
{(xi,Vx,X2,V2,0) e 5R 4 x S" 1 1 = O } on which coordinates 
of attractors (xi, v\, X2, 1*2) are expressed. By employing 
the linear perturbation method which consists of con- 
sidering the solution X(x±, v±, x%, V2) as a superposition 
of a very small perturbation Y{8x\, 5v\, 6x2, SV2) to the 
steady state Vo(a;io, ^io, £20, ^20)) that is X = Vo + Y, we 
obtain the matrix variational equation 

Y = DG(V )Y (8) 

where DG(Vq) is the 4x4 Jacobian matrix 



DG(V ) = 



10 

a -b 

1 

c /3 -0 + c) 



(9) 



describing the vector field along the solution Y and Vo 
being an equilibrium point or steady state with a = 
-l(l+3n(x w -x ));l3 = -'j(l+3Tr(x 2 n-x )). The solu- 
tion of equation (8) after one period T of the oscillations 
in the linearized Poincare map is given by 



Y(T) = Y(0)exp(DG(V )T) 



(10) 



where DG(Vo) represents the time-independent stabil- 
ity matrix of a periodic orbit connecting arbitrary in- 
finitesimal variations in the initial conditions Yq = Y(0) 
with corresponding change Y(T) after one period T. 
The real parts of the roots of the characteristic equation 
det(DG(V ) - IX) = 0, given by 



A 4 + 



At A 3 



AiA + 4o = 



(11) 



determines the stability of the periodic motion, with 
A 3 = 2b + c, A 2 = b(b + c)-(a + f3), A\ = ~b(/3 + a)+ac 
and Aq = a(3. Here A = (Xj) represent the eigenvalues of 
DG(Vq). If one assumes xq rj xiq ~ X20, as in our numer- 
ical experiment, then a = = — In general, complex 



eigenvalues occur in complex conjugate pairs. Let us con- 
sider ctfc, Pk as the real and the imaginary parts of Afc, 
respectively (Afc = ak+iPk)- If Afc is real, the eigenvalues 
are simply the rate of contraction {a% < 0) or expansion 
(afc > 0) near the steady state. If Afc is complex, its real 
part afc gives the rate of contraction (afc < 0) or expan- 
sion (afc > 0) of the spiral while its imaginary part (3k 
contributes for the frequency rotation. The eigenvalues 
of the linearized Poincare map may thus be written as : 

a k = cxp(afcT) (cos(/3 fe T) + % sin(/3 fe T)) (12) 

It turns out from equation (12) that if afc < for all 
Afc , then all sufficiently small perturbations tend forward 
zero as t — > 00 and the steady state (node (n), saddle- 
node (sn), spiral (sp)) is stable. If afc > for all Afc, then 
any small perturbation grows with time and the steady 
state (n, sn, sp) is unstable. In addition if there exist m 
and I such that a m < and a; > 0, the equilibrium state 
is unstable and is called a saddle. Having recourse to the 
above analysis it follows that saddle- node, sn (Afc = +1), 
period-doubling, pd (Afc = — 1), Hopf (f3k 7^ 0, with afc < 
0) and symmetry-breaking (sb) bifurcations are expected 
to occur in the coupled ratchets. Similar bifurcational 
scenario have been found in periodically driven coupled 
single ^l| and double 0] well Duffing oscillators. 
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FIG. 3: Two parameter phase diagram in the (a — uj) plane. 
Regions of period-2 (o), period-3 (•), period-4 (+), chaos 
(black) and non-attraction (white) can be identified. 
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IV. BIFURCATION DIAGRAMS 

In order to investigate the dependence of the system 
on a single control parameter (in this case, the amplitude 
a of the driving function), several bifurcation diagrams 
have been computed, some of which have been chosen to 
illustrate the general structure of the system. Each bi- 
furcation diagram shows the projection of the attractors 
in the Poincare section onto the or the vi coordinates 
versus the control parameter. We employ this technique 
in our numerical exploration using the period T = 2n/to 
of oscillation as the stroboscopic time; with the standard 
4th order Runge Kutta algorithm. 

In what follows, we consider the bifurcational prece- 
dencies of formation of the stable synchronized chaos in 
the synchronization regime i.e in the neighbourhood of 
0.89 < c < 1.05. We begin with a general overview of 
the behaviour of the two coupled systems by displaying 
regions of existence of different periodic and chaotic or- 
bits in a two parameter phase diagram which is plotted 
in (a — lo) plane at fixed values of b = 0.1 and c = 0.5 
(fig- EH)- We obtained the parameter phase diagram using 
the software dynamics [2|| . Different regions of periodic 
orbits of period-2 (o), period-3 (•), period-4 (+) as well 
as chaotic orbits (black region) are clearly visible. Be- 
sides the periodic and the chaotic orbits, there also exist 
non-attracting points (white region) for which the orbits 
diverge. 

In fig. 01 we show bifurcation diagrams for a compar- 
atively small driving frequency of w = 0.3. Different 
bifurcation sequences occur as the driving amplitude a 
is varied. For instance in fig. Eta), there is a resonance 
near a = 0.0559 while a period-1 attractor emerges in 
fig- EJb) from chaos, through a reversed period-doubling 
(pd) and subsequently follows the familiar pd route to 
chaos. For much larger amplitude of the driving force, a 
large period-2 is suddenly created at around a = 1.5612 
from a chaotic region, (see fig. Efc)). This period-2 
attractor is finally annihilated in a crisis event around 
a = 1.5638 leading to a sudden chaotic state (transcient 
chaotic state). After this chaotic state we find bubbles 
of bifurcation dominating periodic windows (further, see 
fig. Efa) as a prototype). When the driving frequency lo 
is increased, saddle-node (sn), symmetry breaking (sb) 
bifurcations, quasiperiodicity as well as period-doubling 
(pd) cascade predominate. For instance, fig. OJa) shows 
the occurence of sn in the vicinity of a = 0.05925 and 
a pd cascade, at around a — 0.06075, leading to chaos 
for lo = 0.4. Moreover sb and sn resulting from a re- 
versed pd are clearly visible in fig. EI (b) for to = 0.5 and 
(c) for uo = 0.75. Similar routes identified in this trans- 
port model have also been found in the coupled Duffing 
model El El. 

To complete our bifurcation study, we proceed to in- 
vestigate large values of the driving frequency (to > lo c — 
0.8) and different ranges of a. Here we find that the ma- 
jority of the earlier observed bifurcation scenario are also 
repeated except that for large a, typically greater than 
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FIG. 4: Bifurcation diagrams for the driving frequency lu = 
0.3, showing (a) a resonance near a = 0.0559, (b) a period-1 
and (c) a large windows of period-2, both sandwiched in the 
chaotic regions. 



the critical value a c = 1.17, periodic attractors dominate. 
Fingerprints of such behaviour can readily be observed 
in fig. ]5j[c) for lo = 0.75, where a period-1 attractor is 
created via a sb bifurcation. 



V. CHARACTERIZING CHAOS 

The best-known characteristic of a chaotic system is 
unpredictability, resulting from sensitivity to initial con- 
ditions (SIC). The slightest deviation from a given set 
of initial conditions yields a completely different output. 
In fig. two typical chaotic trajectories originated from 
two slightly different set of initial conditions (x2 = 0.0 
for A and X2 — 0.01 for B) are displayed. The param- 
eters have been selected from the chaotic region of the 
phase diagram of fig. with the driving amplitude kept 
fixed a = 0.08092844. Although both trajectories ini- 
tially evolve together, after several time steps the state 
variables become essentially independent. This is a clear 
manifestation of the SIC, as witness of chaos. 
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FIG. 5: Bifurcation diagrams for the driving frequencies (a) 
ui — 0.4, (b) uj — 0.5 and (c) ui = 0.75 showing clearly the 
occurence of sn, sb, reversed pd and pd route to chaos. 



Though SIC has come to be the hallmark of chaos, 
it is not unique to the chaotic regime. For example, a 
similar behaviour can be observed in the intermittency 
transition regime where the time series is nearly periodic. 
Thus, to better validate the results obtained above, it is 
necessary to investigate the dynamics of the ratchets in a 
more concrete space - the poincare cross-sections. Anal- 
ysis of the dynamics in such space may help to eventu- 
ally uncover the hidden structure in the behaviour of the 
coupled ratchets. In the poincare cross-section, all tra- 
jectories for a given set of control parameters converge to 
a single object (the attractor) in spite of SIC. The con- 
vergence to an attractor is guaranteed even if the control 
parameters are slightly varied. 

The set of Lyapunov exponents provides an intu- 
itively appealing and yet a very powerful measure of SIC 
and dissipation, both of which are required for a chaotic 
system. All X/. originates from linear stabilty analysis. In 
this approximation, all solutions are of the form exp(Afci), 
k = 1,2, ...M. The maximum Lyapunov exponent X ma x 



km -ln(||L(r)||), 



(13) 



where ||L(r)|| = {5x\ + 8v\ + 8x\ + dv^)? is obtained 
in the Poincare cross-section by solving numerically the 
variational equation (8) simultaneously with the systems 
(4) and (5). A positive Lyapunov exponent is a signa- 
ture of chaos while zero and negative values of the expo- 
nent is an indication of a marginally stable or quasiperi- 
odic orbit and periodic orbit, respectively. In kg. EJb), a 
spectrum of the X m ax, corresponding to the bifurcation 
diagram of kg. \&a) obtained for u> = 0.3, is displayed. 
Regions of ckaotic and periodic solutions are well charac- 
terized. Various bifurcation scenario can be found includ- 
ing sn, sb, quasiperiodic region and pd cascade leading to 
chaos. Interestingly, bubbles of bifurcation beneath the 
chaotic region at a = 1.5630, which by eyes test seems 
to be chaotic as well is not the case indeed. The Lya- 
punov spectrum shows that this regime is just a purely 
quasiperiodic one which yields precisely a period-4 at- 
tractor through a reversed pd. 



1000 2000 3000 4000 5000 6000 7000 

FIG. 6: Typical chaotic trajectories of particle for two slightly 
different sets of initial conditions (a!i = vi = X2 = V2 = 0.0 
for A and x± — vi — V2 = 0.0, X2 = 0.01 for B). This is a clear 
manifestation of SIC. 

Finally in kg. 03 we visualise the attractors in the 
Poincare cross-sections for a kxed valued of tke driving 
force a = 0.080928844 and for several values of the driv- 
ing frequency u> taken from kg. [31 We knd tkat as u> is 
increased, tke ckaotic attractor skrinks (See fig. El (a) f° r 
uj = 0.4, (b) for w = 0.5 and (c) for w = 0.75). As a con- 
sequence, tke attractor size gets enlarged as u> decreases. 
Essentially, this phenomenon has been attributed to tke 
collision of tke attractor witk a periodic orbit in tke ex- 
terior of its basin and is tkus referred to as boundary 
crises 27] . For larger u> values, quasiperiodic and peri- 
odic attractors were found to dominate. As an example, 
figEId) skows a quasiperiodic attractor of period-2 for 

LO = Ul c = 0.8. 
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FIG. 7: (a) Bifurcation diagram for u> = 0.3 and (b) the 
corresponding Lyapunov spectrum in the Poincare map. 
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FIG. 8: Typical attractors for parameters taken from fig. 3, 
with a fixed value of the driving amplitude a — 0.08092844. 
Shown are chaotic attractors (a) for lo — 0.4, (b) for lo — 0.5, 
(c) for lo — 0.75 and (d) a quasiperiodic attractor of period-2 
for lo = 0.8. 



VI. CONCLUSIONS 



In summary, we have investigated the dynamics 
of unidircctionally coupled deterministic ratchets and 
have shown varieties of bifurcation sequences including 
quasiperiodic route to chaos. Using the standard method 
of linear stability analysis, we have examined the stabil- 
ity of the steady state solution of the system leading to 
different types of bifurcations likely to occur in the neig- 
bourhood of the synchronized region. For a given driv- 
ing frequency, the bifurcations depend strongly on the 
values of the driving amplitude a and are generally com- 
plicated. Besides quasiperiodicity, the familiar period- 
doubling and crisis route to chaos were also observed. 
In addition, the coupled ratchets exhibits symmetry- 
breaking bifurcations, resonance, saddle-node and bub- 
bles of bifurcations. Using Lyapunov exponent spectrum 
as well as sensitivity to initial conditions, we character- 
ized chaos in this system. A perusal of the Poincare 
cross-sections revealed boundary crises in which the size 
of chaotic attractor is suddenly enlarged as the driving 
frequency is gradually decreased for a fixed driving am- 
plitude. Although this system presents a very chaotic 
structure, it remains nevertheless ordered as the driving 
frequency becomes large, thereby providing critical pa- 
rameters for a chaotic transport of particles. This makes 
the present study essentially interesting for technological 
applications. 
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